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1  Introduction 

There  has  been  significant  interest  in  recent  years  in  the  combustion  of  sub-micron  aluminum 
particles,  and  the  standard  theoretical  framework  is  the  shrinking-core  model.  This  model,  in  the 
context  of  any  fuel,  has  roots  that  go  back  60  years,  [1].  A  common  reference  is  a  chapter  in 
[2],  although  the  discussion  there,  within  a  framework  of  global  statements,  is  reminiscent  of  well- 
stirred  reactor  discussions,  familiar  in  the  chemical  engineering  literature,  so  that  certain  details 
are  obscured.  However,  the  model  is  simply  described:  A  spherical  core  of  liquid  aluminum  is 
surrounded  by  a  shell  of  alumina  and  the  aluminum  core  shrinks  as  aluminum  is  converted  to 
alumina  at  the  metal/oxide  interface.  The  conversion  occurs  because  of  the  inward  transport  of  O 
atoms  (we  assume  that  the  surrounding  atmosphere  is  air)  from  the  outer  boundary  of  the  oxide 
shell.  It  is  commonly  assumed  that  this  transport  is  solely  diffusive  in  nature,  and  quasi-steady,  at 
best  an  approximation;  and  that  apart  from  the  aluminum/oxide  transformation,  the  geometry  is 
fixed.  Provided  these  approximations  only  lead  to  modest  error  it  might  be  reasonably  argued  that 
it  is  acceptable  in  view  of  two  significant  uncertainties:  the  value  of  the  diffusion  coefficient  of  the 
O  atoms;  and  the  value  of  the  O  concentration  within  the  oxide  at  the  oxide  boundary.  However,  as 
we  shall  see,  one  of  the  neglected  ingredients  has  significant  mechanical  consequences  in  the  context 
of  the  spherical  geometry,  and  therefore  it  is  important  to  consider  it. 

Because  of  the  mechanical  consequences,  an  analytical  treatment  in  the  spherical  geometry  is  not 
possible,  and  yet  analytical  treatments  are  of  great  value  in  revealing  fundamental  physics.  And  so 
we  shall  start  our  discussion  with  a  planar  model,  purely  because  of  the  insights  thereby  achieved, 
not  because  it  corresponds  to  a  configuration  to  be  found  in  the  physical  world. 

2  The  Planar  Problem 

The  configuration  is  sketched  in  Figure  1.  Liquid  aluminum  occupies  the  interval  0  <  x  <  xai ., 
alumina  occupies  the  region  xai  <  x  <  xoxd.-  The  interface  between  the  two  materials,  also  the 
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reaction  front,  moves  to  the  left  with  speed  vrf.  the  oxide  moves  as  a  solid  body  to  the  right  with 
speed  v.  The  M  are  mass  fluxes  (of  O  atoms,  aluminum,  oxide)  to  or  from  the  front,  measured 
relative  to  the  front. 

This  model  does  not  account  for  aluminum  diffusion,  only  diffusion  of  the  O  atoms.  If  diffusion 
of  both  occurs  the  reaction  would  occur  interior  to  the  alumina  either  at  a  front  or  within  a  spread- 
out  reaction  zone.  There  are  various  references  to  (and  opinions  about)  aluminum  diffusion  in  the 
literature. 

Campbell  et  al,  [3]  carry  out  molecular  dynamics  simulation  in  which  they  consider  a  ball 
of  aluminum  atoms  of  diameter  20nm  exposed  to  an  oxygen  atmosphere  and  observe  oxide  layer 
growth  with  diffusion  of  both  aluminum  and  oxygen.  They  claim  that  the  aluminium  diffusion  is 
greater  than  that  of  the  oxygen  because  of  the  smaller  steric  size  of  the  aluminum  atoms.  But 
the  temperature  is  300K  so  that  the  aluminum  is  solid,  and  here  we  are  concerned  with  liquid 
aluminum;  in  the  rocket  business  the  interest  is  in  aluminum  particles  released  into  the  hot  rocket 
chamber. 

Rai  et  al,  [4]  claim  that  when  the  aluminum  is  solid,  oxygen  diffusion  alone  occurs;  when  the 
aluminium  is  liquid  both  aluminum  and  oxygen  diffusion  occur.  For  evidence  they  note  that  they 
get  hollow  post-burn  oxide  particles  when  the  temperatures  are  high  enough  for  the  aluminum  to 
be  in  the  molten  state.  This  is  an  interesting  argument,  but  if  one  carries  out  a  standard  shrinking- 
core  model  analysis,  with  an  interior  reacting  front,  and  invokes  Ockham’s  Razor  -  the  volume 
of  the  aluminum  sphere  changes  precisely  to  accommodate  the  outward  aluminum  diffusion  -  this 
leads  naturally  to  a  solution  sans  interior  void.  To  get  a  void  one  would  need  to  move  the  outer 
boundary  of  the  aluminum  inward  more  slowly  but  how  one  would  specify  that  modification  is 
unknown.  Note  that  an  appeal  to  the  steric  argument  to  justify  that  alumimim  diffusion  will  occur 
ignores  other  physical  ingredients.  The  atoms  at  a  liquid  surface  are  tied  by  atomic  forces,  and 
an  atom  will  only  detach  from  the  surface  when  the  sum  of  the  collision  forces  from  its  neighbors 
is  sufficient;  sufficiency  is  a  function  of  the  temperature.  And  should  a  central  void  be  inclined  to 
arise  the  pressure  there  will  be  small,  provided  by  aluminum  vapor,  and  there  will  be  a  positive 
radial  pressure  gradient  in  the  liquid  which  will  also  discourage  surface  detachment  of  atoms. 

Trunov  et  al,  [5]  divide  the  oxidation  process  into  two  stages.  At  early  times  the  oxide  is 
amorphous  in  nature  and  is  generated  by  aluminum  diffusion.  But  the  amorphous  layer  is  only 
stable  for  thickness  up  to  5nm,  and  after  that  the  7  crystalline  form  is  generated  (crystal  size  5nm) , 
generated  by  O  diffusion.  The  a  crystal  state  might  occur  later  (crystal  size  45nm)  but  only  for 
large  drops.  However,  these  conclusions  are  rooted  in  experiments  for  slow  oxidation,  not  the  fast 
oxidation  that  is  of  interest  to  us. 

Our  conclusion  is  that  the  role  of  aluminum  diffusion  for  our  problem  is,  as  yet,  unknown,  and 
our  choice  at  this  time  is  to  ignore  it.  In  a  later  study  we  shall  revisit  this  question,  partly  because 
of  the  possibility  that  aluminum  diffusion  plays  a  role  in  transition  solutions  leading  to  the  classical 
drop  burning  regime. 

Returning  to  Figure  1,  the  two  speeds  are  measured  in  the  laboratory  frame  tied  to  x  =  0.  The 
speed  v  arises  from  mass  conservation  at  the  reaction  front  where,  overall, 

2AI  +  30  — y  AI2O3 
(54)  (48)  (102) 

(We  shall  henceforth  use  a=54,  /3=48,  and  7=102  so  that  our  results  can  be  understood  in  other 
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Figure  1:  Configuration,  v  and  vrf  are  measured  in  the  laboratory  frame;  all  M  are  measured 
relative  to  the  reaction  front. 
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contexts.)  Then  if  in  a  frame  attached  to  the  front  the  mass  flux  of  0  to  the  front  is  M0,  the  mass 
flux  of  aluminum  to  the  front  is  Mai.,  and  the  mass  flux  of  oxide  away  from  the  front  is  Moxd.  (then 
all  the  M  are  positive)  we  have 


Mai.  a  , 

=  —  and 


Mn 


p 


Moxd. 
Mai . 


1 

a 


Now 


Mal .  PalVrf 


(2) 


(3) 


and 


M0xd.  —  iy  T  Vrf'jPoxd.  (4) 

Also,  accounting  for  both  convection  and  diffusion  of  0, 

dc 

M0  =  —(v  +  vrf)c  +  D—  (5) 

where  c  is  the  concentration  of  0  atoms  (in  dimensions  that  of  density).  Here  we  shall  assume  that 
the  reaction  is  fast,  and  is  controlled  by  diffusion,  so  that  at  the  front  c  =  0.  Then 


dc 

Mo  |  fast  chem.  =  ^7^ 


Within  the  oxide,  c  satifies  the  equation 


dc  dc  d2c 

dt  V  dx  dx2 


(6) 


(7) 


Note  that  we  neglect  convective  fluxes  of  c  relative  to  the  oxide  which  could  only  arise  because  of 
pressure  gradients.  Reference  4  accounts  for  pressure  gradients,  motivated  by  the  simulations  of 
Reference  3,  and  we  shall  say  something  about  general  stresses  in  a  later  section  when  we  consider 
the  spherical  geometry. 

Insofar  as  the  oxide  speed  v  is  concerned,  we  have  yet  to  come  across  any  explicit  discussion 
of  it,  although  its  physical  origins  and  one  of  the  important  roles  that  it  plays  have  long  been 
recognized  by  corrosion  scientists.  If  a  thin  layer  of  metal  is  converted  to  an  oxide  layer  of  greater 
volume,  as  in  the  case  of  aluminum,  this  converted  layer  will  push  out  previously  generated  layers. 
And  in  the  case  of  curved  displaced  layers,  hoop  stresses  will  be  generated.  Corrosion  scientists  are 
only  interested  in  thin  surface  layers,  [6,  7,  8].  More  recently  this  phenomenon  has  been  recognized 
by  those  concerned  with  aluminum  drop  ignition,  [9],  and  here  also  it  is  only  thin  layers  of  surface 
oxide  that  are  of  interest.  Our  interests  go  deeper  than  that,  literally. 

It  should  be  noted  that  Ref  [9]  proposes  a  model  of  sub-micron  aluminum  combustion  that  is 
quite  different  from  the  shrinking-core  model.  It  is  rooted  in  the  idea  that  genesis  layers  of  oxide 
violently  rupture  when  heated  rapidly,  leading  to  an  expansion  wave  which  propagates  into  the  core 
so  that  liquid  aluminum  is  violently  expelled  in  atomic  sized  clusters  of  naked  aluminum  which  is 
rapidly  consumed.  The  fundamental  argument  is  a  thermal  one,  quite  independent  of  combustion, 
and  Lynch  et  al,  [10],  subjected  particles  to  the  high  temperature  behind  a  reflected  shock  in  an 
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inert  (argon)  atmosphere  and  spectrally  examined  the  post  shock  gases.  They  found  no  trace  of 
aluminum  vapor  for  temperatures  below  2300K. 

Returning  to  equation  (7),  for  the  moment  we  shall  adopt  the  quasi-steady  approximation, 
neglecting  the  time  derivative.  Provided  dynamical  ingredients  are  not  at  the  heart  of  a  physical 
problem  this  can  often  be  a  useful  approximation,  providing  valuable  insights,  albeit  with  quan¬ 
titative  error.  We  shall  retain  the  time  derivatives  in  numerical  calculations  reported  later  in  the 
paper. 

The  general  solution  is 


c  =  A  +  B  exp 


for  arbitrary  constants  A  and  B  to  be  determined  from  the  boundary  conditions 


(8) 


c  =  0  at  x  =  xai  and  c  =  c0  at  x  =  xoxd.  (9) 

As  noted  earlier,  the  determination  of  cG  is  difficult.  O  atoms  are  generated  by  collision  of  oxygen 
molecules  with  the  aluminum  surface,  collide  themselves  with  the  surface,  and  a  certain  fraction 
attach  and  subsequently  pass  into  the  alumina  between  the  alumina  molecules.  The  gas  phase 
processes  are  not  in  the  continuum  domain  for  there  is  a  Knudsen  layer  at  the  drop/particle  surface. 
With  c  then  specified  as  a  function  of  v,  D,c0,xai.  and  xoxd.  equations  (6),  (3)  and  (2)  yield 


Now 


M0  =  vc0 


exp  -  (aw.  -  xai.)  -  1 


P 

Pal.Vrf 

a 


where,  from  (4),  (3),  and  (2) 


dxgi. 

dt 


—vrf  and 


dxoxd. 


dt 


=  V 


(10) 


(11) 


V  = 


1  Pal. 
®  Poxd. 


Vrf 


(12) 


With  pai= 2400  kg/m3,  /0Oa;d.=3950  kg/m3  the  constant  coefficient  on  the  rhs  is  0.148  for  the  alu¬ 
minum  problem. 

Equations  (12),  (11),  and  (10)  define  equations  for  xai  and  xoxd.  from  which  we  can  deduce 


where 


Xoxd.  Xoxd.  |o  —  [Xai.  ^aijo] 


7  Pal. 
&  Poxd. 


(13) 


%oxd.  |o  ^al.\o  ^ 

the  initial  thickness  of  the  aluminum.  And 


(xaL  -  L)2  =  2A t 


(14) 


(15) 


where 
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A  =  Dr-^l  +  T]"1!!! 


(16) 


r 


i  +  r 

P  Pal. 


7  Pal. 
&  Poxd. 


Since  ca  is  three  orders  of  magnitude  smaller  than  pai .,  the  logarithm  function  can  be  replaced  by 
a  linear  approximation,  whence 


A~  D(1  +  r)"1^—  (17) 

P  Pal. 

At  the  end  of  the  burn,  when  xai=  0,  t  =  tburn  oc  L2  a  ”d2  —  t”  law  in  the  vocabulary  of  fuel-drop 
burning.  The  formula  is 


tb 


urn 


L2  15  Pal.  (1  +  r) 
2  a  c0  D 


(fast  chemistry,  quasi-steady) 


(18) 


A  d2  —  t  law  arises  for  a  simple  dimensional  reason  for  the  only  way  to  define  a  time  is  L2  /  D  (cf. 
(15),  (16).  As  noted  earlier,  the  specification  of  D  is  challenging,  but  if  the  burn-law  is  essentially 
a  d2  —  t  law  an  order-of-magnitude  estimate  for  D  is  T2/tburn-  an  estimate  that  will  be  valid  for 
the  spherical  geometry  with  L  equal  to  the  particle/drop  diameter. 

For  the  same  dimensional  reasons  the  d2—t  law  arises  in  the  classical  model  of  external  flame  fuel- 
drop  burning  in  which  the  flame  is  essentially  described  by  the  Burke-Schumann  limit.  Sufficiently 
large  aluminum  drops  will  essentially  burn  in  this  fashion,  modified  slightly  by  the  presence  of 
an  alumina  cap  on  the  aluminum  ball.  A  thorough  review  of  super-micron  sized  drop  burning  is 
presented  in  [11],  and  Figure  1  of  that  paper  shows  burning  times  vs  diameter  from  12  experimental 
studies,  with  diameters  ranging  from  about  16  microns  to  800  microns.  There  is  substantial  scatter 
in  the  combined  data,  but  a  best  fit  corresponds  to  a  d1"  —  t  law.  We  know  that  when  the  flame 
temperature  is  greater  than  the  ambient,  the  theoretical  burning  response  of  maximum  temperature 
vs  Damkohler  number  (proportional  to  d2)  is  S-shaped  with  most  observed  burning  occurring  on  the 
upper  branch.  Then  a  reduction  of  the  Damkohler  number  will  lead  to  quenching  when  the  upper 
branch  ends  unless  the  ambient  temperature  is  large  enough  to  sustain  significant  evaporation  of  the 
aluminum  so  that  combustion  can  continue  on  the  lower  branch.  Thus  at  an  ambient  temperature  of 
1500K,  for  example,  sufficiently  small  drops  can  not  burn  classically;  on  the  other  hand  submicron 
drops  can  burn  at  that  temperature  in  an  internal  combustion  context.  The  nature  of  the  transition 
zone,  whether  or  not  one  always  exists  even,  is  not  understood  at  the  present  time.  Lynch  et  al, 
[12],  examined  3-11  micron  drops  at  temperature  greater  than  2400K,  and  therefore  greater  than 
the  melting  temperature  of  alumina.  They  used  a  variety  of  atmospheres  and  obtained  dn  —  t  laws 
for  n  between  1.3  and  0.2. 

Returning  to  the  solution  summarized  by  equations  (15)  and  (16)  we  can  write  down  formulas 
for  the  speeds.  We  have 


and 


(19) 
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V  = 


(20) 


1  Pal. 


-  1 

®  Poxd. 

The  singularity  at  t  =  0  arises  because  at  t  =  0  we  have  c  =  cQ  at  x  =  xoxd.\o  and  c  =  0  at 
x  =  xai, |o,  coincident  locations.  It  is  the  familiar  singularity  associated  with  the  Rayleigh  problem, 
for  example. 

If  equation  (19)  is  integrated  over  the  burn  time,  the  final  location  of  the  outer  oxide  boundary 
is  easily  calculated 

|  _  T  Pal.  r  foi  ^ 

Xoxd.\burn 

Ot  Poxd. 

This  simply  states  that  the  final  mass  of  oxidizer  poxd.xoxd.  \bum  is  equal  to  ^  times  the  initial  mass 
of  aluminum  pai.L,  as  required  by  the  stoichiometry  of  (1).  If  the  convection  term  in  the  equation 
for  c,  (7),  is  neglected,  such  global  conditions,  at  different  times,  can  bypass  any  discussion  of  v. 
Global  mass  considerations  can  also  lead  to  the  prediction  of  the  final  drop/particle  radius  for  the 
spherical  problem,  but  there  the  role  of  v  is  more  profound  because  the  generation  of  a  radial  oxide 
speed  at  the  aluminum/alumina  interface  can  not  generate  a  mere  solid  body  displacement,  but 
must  lead  to  a  significant  stress  field.  This  is  the  subject  of  a  later  section. 

Consider  now  the  importance  of  the  convection  term  in  (7).  Since  x  ~  L,  c  ~  c0,  and  v  v/D Jt 

we  have 

dc  c"  '  dS-5:  (22) 


dx  t  “d  t 

and  convection  may  not  be  neglected.  But  the  time  derivative  has  the  same  estimate,  and  so  should 
be  retained  also;  done  properly,  this  is  a  2-boundary  Stefan  problem  with  convection.  Fortunately 
dynamics  is  not  central  to  the  problem,  and  this  would  not  be  the  first  time  that  a  quasi-steady 
model  has  been  adopted  in  circumstances  for  which  it  is  but  an  extra-rational  approximation. 
We  can  be  reasonably  sure  that  the  qualitative  picture  that  we  have  painted  here  is  essentially 
sound,  and  the  omitted  time  derivatives  provide  but  quantitative  corrections.  Some  measure  of 
the  magnitude  of  these  corrections  will  be  provided  by  numerical  calculations  described  in  the  next 
section. 

We  finish  this  section  by  identifying  the  burn  time  when  v  is  neglected,  and  xoxd.  remains  fixed 
at  the  initial  value  L.  Then  the  solution  for  c  is  linear  in  x,  vanishing  at  xai  and  having  the  value 
c0  at  L.  Thus  the  diffusive  flux  to  the  reacting  front  is  a  simple  function  of  xai,  proportional  to 
vrf.  In  this  way  we  find  that  the  burn  time  is 


^bum  =  —  —  (fast  chemistry,  quasi-steady,  v  =  0) 

2  a  D  cQ 


(23) 


which  differs  from  (18)  by  the  factor  (1  +  T).  In  section  4  we  compare  this  burn  time  and  the  burn 
time  defined  by  (18)  with  numerical  results  obtained  with  v  and  dc/dt  retained. 


3  Finite  Chemistry  Effects 

In  this  section  we  start  with  the  model  of  section  2  but  modify  it  to  account  for  finite-rate  chemistry 
at  the  reaction  front,  so  that  c  does  not  vanish  there.  The  aluminum  is  liquid  but  how  one  should 
characterize  the  ensemble  of  O  atoms  is  not  clear.  The  constraints  of  the  alumima  molecules  rob 
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it  of  the  essential  characteristics  of  a  gas;  nor  is  it  a  liquid.  And  so  the  reaction  should  probably 
be  thought  of  as  of  the  heterogeneous  kind. 

Global  modeling  of  homogeneous  gas-phase  reactions  has  proven  to  be  most  valuable.  If  there 
are  two  reactants  concentrations  A,  B,  which,  in  a  global  sense  react  according  to 

aA  +  /3B  — >  Products  (24) 

(i.e.  this  is  not  an  elementary  reaction)  then,  with  a  nod  to  the  law  of  mass  action  and  to  the 
Maxwellian  distribution  of  molecular  speed  (related  to  collision  energy)  we  can  model  the  reaction 
rate  by 


K  =  kA^B1'  exp 


E 

rT 


(25) 


with  suitably  chosen  values  of  k,  and  E. 

This  can  also  work  well  for  certain  heterogeneous  reactions;  for  example,  experimental  data  for 
the  chemical  erosion  of  carbon  by  hot  gases  (relevant  in  nozzle  erosion)  has  been  fitted  with  similar 
formulas:  a  discussion  with  several  references  may  be  found  in  [13].  This  is  not  too  surprising. 
In  one  case  gas  molecules  are  colliding  with  gas  molecules  with  a  collision  frequency  that  depends 
on  the  two  concentrations;  in  the  other,  gas  molecules  are  colliding  with  the  surface  of  the  carbon 
with  a  collision  frequency  that  depends  on  the  single  concentration.  In  both  cases  the  energy  of 
the  collisions  is  related  to  the  high  temperature  tail  of  the  Maxwellian. 

We  know  of  no  authoritative  discussion  of  a  global  kinetics  model  for  our  problem.  However,  it 
is  not  implausible  to  speculate  that  the  rate  at  which  aluminum  atoms  encounter  O  atoms  depends 
on  the  concentration  of  the  latter  and,  since  we  do  not  account  for  thermal  effects  so  that  the 
temperature  dependence  of  the  efficacy  of  such  encounters  is  immaterial,  we  can  take 


K  =  kcv  at  x  =  xai  (26) 

to  be  the  reaction  rate.  That  is,  with  the  choice  rj  =  1,  equation  (6)  is  rewritten 

b  dc 

M0  ~ Pal.^rf  {v  T  Vrf)crf  T  D  |  rf  kcrj  (27) 

With  this  modification,  the  analysis  proceeds  as  in  Section  2,  so  that  we  do  not  need  to  burden  the 
reader  with  the  details.  The  results  are: 


where,  as  before, 


f'LU.  ' 

OL  OL  Poxd. 


(V  \ 

(v  W 

exp \DXaL) 

-  exp  xoxd.j 

V  =  V. 


f  / 


7  Pal. 


-  1 


&  Poxd. 


vrf  — 


dxoxd. 

dt 

dxal . 


(28) 


(29) 


dt 


To  examine  equation  (28)  it  is  useful  to  introduce  some  scalings.  To  this  end  we  write 


D 


D 


•Eoxd.  LSoxd.i  ■Kal.  V  j  ^ rf 


so  that  (28)  becomes 


(30) 


Co 


13  Pal. 

D 

1 

7  Pal. 

a  c0  kL 

OL  poxd. 

'p  ^ 

7  Pal. 

—  Pal .  + 

- C 

Ya 

Cl  poxd. 

exp  (usoxd.)  - 


O  [exp(usaZ.) 


exp  (usaL)  = 
-  exp(usra(i)] 


(31) 


In  the  fast  chemistry  limit  we  recover  the  results  of  section  2.  In  the  slow  chemistry  limit  (k  — >•  0) 
u  vanishes  so  that  the  exponentials  can  be  simplified,  and  (31)  leads  to 


Urf 


PR 

a  kL 


T  ( Soxd .  ®aZ.) 


P  +  7  Pal.  Cp 
OL  a  poxd.  Pal. 


which,  when  D/kL  1,  leads  to 


Co 

Pal. 


(32) 


Urf  ~ 


a  c0  kL 
13  Pal.  D 


i.e.  vrf  ~ 


OL  Cq 

« — k 
P  Pal. 


(33) 


Then  crf  =  ca  (there  is  no  gradient  of  oxygen  in  the  oxide),  and  the  burning  law  is  a  linear  one; 
the  only  relevant  speed  is  k: 


tb 


urn 


L 13  Pal. 
acQk 


(slow  chemistry) 


(34) 


This  familiar  result  mirrors  Glassman’s  discussion  for  gas- phase  burning,  [14].  Note  that  quite 
independently  of  the  simple  limit  solution  for  c,  the  scalings  show  that  the  ratio  of  the  diffusion 
term  to  both  the  time  derivative  and  the  convection  term  is 


O 


D  Pal. 
kL  c0 


and  so  the  neighborhood  of  this  limit  may  be  discussed  within  a  purely  diffusive  framework. 


4  Numerical  Results 

In  section  2  we  deduced  that  the  quasi-steady  burn  time  in  the  fast  chemistry  limit  is  L2/ 2A 
(equation  (18));  in  the  same  limit  but  with  v  neglected  it  is  L2 f3pai./2aDc0  (23);  and  in  section  3 
we  deduced  that  the  burn  time  in  the  slow  chemistry  limit  is  (/ 3Lpai  )/(akc0 )  (34).  Consider  the 
parameter  choices 


c0  =  6.4kg/m3  k  =  18.75m/s  D  =  10  5m2/s  pai  =  2400kg/m3  poxd.  =  3950kg/m3 

a  =  54  (3  =  48  7  =  102 


(35) 
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Then  the  limit  burn  times  are 


1.91  x  107L2  s  eq(18)  1.67  x  107L2  s  eq(23)  17.78L  s  eq(34)  (36) 

In  the  log-log  plane  the  first  and  third  of  these  define  two  straight  lines  which  intersect  at  L  = 
.931  x  10-6,  with  kL/D  =  1.75;  order  1  values  of  kL/ D  characterize  the  transition  from  the 
linear  burn  law  to  the  quadratic  law.  This  transition  is  more  accurately  described  using  numerical 
simulations  in  which  the  time  derivatives  are  retained. 

In  doing  this  we  note  that  the  system  is  stiff,  because  of  the  small  value  of  c0/pai .,  so  that  an 
explicit  solver,  such  as  Runge-Kutta,  is  inadequate.  We  used  a  Crank-Nicholson  scheme,  and  the 
results  are  shown  in  Figure  4.  The  straight  lines  cross  at  L  =  10~6,  kL/D  =  1.9,  slightly  different 
from  the  conclusions  using  the  analytical  results.  At  L  =  ICR9  the  analytically  derived  burn  time 
is  1.78  x  10”8  in  close  agreement  with  the  numerical  result.  At  a  burn  time  of  1  the  analytically 
derived  value  of  L  from  eq(18)  is  2.29  x  10  ,  from  eq(23)  it  is  2.45  x  10-4.  In  both  cases  the  errors 

are  quite  small. 

5  Summary 

In  the  past  sections  we  have  examined  the  familiar  shrinking-core  model  in  the  planar  context  in 
order  to  highlight  several  features  that  do  not  appear  to  have  been  adequately  discussed  in  the 
past.  One  is  the  speed  v  acquired  by  the  alumina  because  of  the  addition  of  mass  (aluminum  and 
O)  at  the  reaction  front.  It  should  be  noted  that  if  Rai  et  al,  [4],  are  correct  in  arguing  that  the 
reacting  front  or  domain  lies  within  the  alumina,  then  there  will  be  convective  fluxes  of  alumina 
away  from  the  domain  on  both  sides  as  alumina  is  added  by  reaction,  necessarily  tearing  apart 
existing  alumina  molecules,  and  there  would  be  substantial  material  stresses  associated  with  this. 
But,  as  we  shall  see  in  the  next  section,  there  are  substantial  stresses  associated  with  the  creation 
of  alumina  in  the  context  of  the  model  that  we  adopt  when  the  geometry  is  spherical  so  that  simple 
solid-body  displacement  of  the  alumina,  as  in  the  planar  problem,  is  not  possible. 

A  second  point  that  we  have  discussed  is  the  role  of  the  convection  term  and  the  time  derivative 
in  defining  the  distribution  of  c,  the  O  concentration.  Technically,  both  should  be  retained  in  the 
fast-chemistry  limit,  as  both  terms  are  comparable  to  the  diffusion  term.  However,  the  time  scale 
is  defined  by  vrf,  not  v,  and  for  aluminum  burning  v  is  roughly  1/7  of  vrf.  Then  an  error  of  15% 
caused  by  omitting  the  convection  term  may  not  be  of  much  significance  if  D  is  not  well  known. 
There  is  no  good  a  priori  reason  for  omitting  dc/dt,  but  there  also  the  error  is  small,  at  least  for 
the  values  used  in  aluminum  combustion. 

6  The  Spherical  Problem 

In  this  section  we  turn  from  the  planar  geometry  to  the  spherical  one.  Here  equations  (l)-(6)  are 
still  applicable,  as  is  (7)  when  the  diffusion  term  is  written  in  spherical  form,  but  the  oxide  speed  v 
presents  difficulties,  since  it  no  longer  corresponds  to  a  solid-body  displacement  unless  we  suppose 
that  the  alumina  is  rigid.  To  calculate  v  we  have  to  consider  the  displacement  and  stress  fields  in 
the  oxide,  and  that  is  our  focus. 

Consider  a  sphere  of  liquid  aluminum  surrounded  by  a  shell  of  alumina.  As  before  we  assume 
that  the  reaction  surface  is  at  the  boundary  between  the  two  materials  without  assumption  as 
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Figure  2:  Burn  Times 
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to  what  controls  the  reaction,  whether  it  be  diffusion  of  O  through  the  oxide,  or  reaction  at  the 
aluminum  surface.  We  assume  that  the  reacting  surface  moves  inwards  with  speed  vrf,  and  that 
the  creation  of  alumina  at  the  surface  leads  to  an  outward  flux  of  alumina  from  the  surface  with 
local  speed  v.  We  calculate  the  stress  consequences  of  this  by  considering  an  elastic  shell  of  outer 
radius  b  and  inner  radius  a.  a  decreases  with  time  because  of  the  aluminum  consumption,  and  we 
shall  assign  this  function,  b  changes  for  two  reasons:  the  volume  of  oxide  generated  is  greater  than 
the  volume  of  aluminum  consumed;  elastic  displacement. 

Considering  just  the  first,  if  the  initial  radii  are  bo  and  ao,  and  the  final  value  of  b  is  bf  (a/  =  0) 
then 

b'  =  b°+{^)a°-  (37) 

bo  is  greater  than  ao  by  an  amount  s,  the  thickness  (small)  of  the  genesis  oxide  layer.  Then  it  is 
trivially  shown  that 

/  <  1.047  (38) 

bo 

and  we  shall  neglect  this  change. 

As  for  the  elastic  displacement,  our  calculations  are  within  the  framework  of  linear  elasticity, 
and  so  that  is  also  neglected. 

The  pressure  in  the  external  environment  is  taken  to  be  Pb  and  that  in  the  core  region  Pa.  The 
presence  of  aluminum  in  the  core  plays  no  role  except  that  its  contribution  to  the  reaction  provides 
a  relation  between  vrf  and  v,  as  before.  In  what  follows  r  is  the  radius,  a  is  Poisson’s  Ratio,  E  is 
Young’s  Modulus,  Tr  is  the  radial  stress,  t\  and  72  are  hoop  stresses,  5  is  the  radial  displacement, 
ei  and  £2  are  the  hoop  strains,  and  er  is  the  radial  strain. 

The  solutions  of  the  linear  equations  of  elasticity  for  this  configuration  are 


1 


Tr  = 


Tl  =  r2  = 


( b 3  —  a3) 

1 

( b 3  —  a3) 


Paa  -  PbbA  + 


Pn.as  -  Puli'  - 


ii(Pb  -  Pa. 

r3 

3fr3(A  -  Pa 
2  r3 


5  = 


E(b3  -  a 3) 


(1  -  2 a)(Paa3  -  Phb 3)  -  l-±^0?b\Ph  -  Pa) 

5  d5 

—  €2  —  er  —  — 

r  or 


(39) 


These  formulas  can  be  derived  in  a  straightforward  fashion  from  the  equations  of  linear  elasticity,  but 
may  be  found  in  [15].  They  were  used  by  Rosenband,  [16],  to  examine  the  stresses  on  the  genesis 
layer  of  amorphous  alumina  created  when  the  aluminum  particle  is  ignited  and  the  aluminum 
expands. 

The  radial  velocity  associated  with  the  displacement  5  is 


d5 

dt 


which  is  equal  to  v  when  evaluated  at  r  =  a 


(40) 


These  solutions  would  normally  be  used  to  calculate  the  displacement  and  stress  fields  when  Pa 
and  Pb  are  assigned.  But  in  our  problem  we  are  interested  in  the  assignment  of  the  radial  velocity 
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at  the  inner  radius  a  and  the  assignment  of  Pb,  with  Pa  to  be  determined  as  the  driver  for  v.  Pa, 
like  a,  will  be  a  function  of  time. 

Rather  than  differentiate  the  formula  for  <5  with  respect  to  t  we  shall  use  a  as  a  surrogate  for 
time,  with 


since 


a(0)  =  b 


a  =  b  —  vrft 

where  we  shall  take  vrf  (and  therefore  v,  see  (12)  )  to  be  constant. 

So  we  differentiate  5  with  respect  to  a  holding  r  constant,  and  then  set  r  =  a.  Thus 


(41) 


Af  +  B4- 

da 


C  = 


E(l-a)  88 
Pb  da 


E(  1  —  a)  f  v  \ 

Pb  \  Vrf J 


where 


(42) 


Pa 

Pb 


a  = 


b 3 


and  the  coefficients  are  given  by  the  formulas 


(43) 


A  =  3a 


«( 1  -  2ff)  +  ^  +  fJ) 


B  = 


3a 


(1  ~a) 


a(  1  -  2cr)  +  -(1  +  a) 


+  3cx(l  —  2cr )  +  — (1  +  <t) 


^  ck  9.,  ,  3,,  . 

c  =  74 - +  c7-) 

(1  —  a)  2  2 


(44) 


Before  we  can  solve  equation  (42)  we  need  to  identify  the  initial  value  of  Pa.  For  this  purpose 
we  shall  suppose  that  there  is  a  genesis  layer  of  amorphous  oxide  of  thickness  5nm  and  that  the 
particle,  of  temperature  300K,  is  placed  in  an  atmosphere  of  temperature  greater  than  or  equal  to 
the  aluminum  melting  temperature.  Because  the  particle  is  small  and  its  thermal  conductivity  high, 
it  is  assumed  that  its  temperature  is  raised  instantly  to  the  atmospheric  value,  and  so  experiences 
a  significant  increase  in  volume  which  imposes  stresses  on  the  alumina.  This  is  a  problem  that  has 
been  examined  before,  [9]. 

For  this  problem  the  thin-shell  approximations  to  equations  (39)  are  valid.  Thus,  writing 


we  have,  for  s  small 


(b  —  a)  =  s,  r  =  a  +  7s 


(45) 


Tr  ~  ~Pa  +  7 {Pa  ~  Pb) 

Tl  =  T2  ~  j.  j  (Pb  ~  Pa) 


a(  1  —  a)  a 


(46) 
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Then  with  E  =  300GPa,  <7=0.25,  S/a=. 04  (corresponding  to  a  liquid  aluminum  density  of  2400kg/m3, 
a  solid  aluminum  density  of  2700kg/3),  and  s/a  =  1/16  we  have 

(Pa  —  Pb)  ~  2GPa,  t\  =  T2  ~  — 5.3GPa  (47) 

Then,  to  solve  (42)  we  take 

—  =  0.148  E  =  300  GPa  a  =  0.25  Ph  =  2  x  106  Pa  Pa  =  2GPa  (48) 

vrf 

where  the  last  two  are  initial  conditions  applied  at  a  =  0.834. 

Figure  6  shows  the  variations  of  cp  with  a/b.  From  this  we  see  that  the  internal  pressure  rises 
significantly  as  the  ball  of  aluminum  shrinks.  In  this  connection  we  note  that  Campbell  et.al.[ 3] 
conducted  a  molecular  dynamics  study  of  the  early  stages  of  the  oxidation  of  solid  aluminum  and 
found  that  significant  internal  pressures  were  generated.  And  Rai  et.al.  [4]  were  motivated  by 
this  work  to  impose  pressure  gradients  in  their  study  of  combustion  with  both  O  diffusion  and  A1 
diffusion. 

Once  cp  is  known  we  can  calculate  the  stresses  and  the  strains.  Figure  6  shows  the  strains  (vs 
a/b )  at  r  =  a  and  r  =  b 


-I 

a 


1  Pb 

(1  —  a)  E 

1  Pb 


l,  (1  -at)E 


(1  -  2a)(cpa  -  1)  -  ^y^(l  -  <P) 
(1  —  2a)(cpa  —  1) - — — a(l  —  cp) 


(49) 


And  figure  6  shows  the  hoop  stresses  at  r  =  a  and  r  =  b. 


T\.2\r=a  — 


T" 1,2  r=b  ~ 


Pb 


(1  -  a) 

Pb 


(1  -  a) 


(j)a  -  1  - -(1  -  (/>) 
(pa  —  1  —  *a(  1  -  <p) 


(50) 


Clearly,  when  a  becomes  sufficiently  small  the  strains  and  stresses  can  not  be  sustained  by 
the  alumina,  and  fracture  must  occur,  albeit  not  in  the  neighborhood  of  the  outer  surface.  The 
breaking  of  spherical  symmetry  in  this  way  must  have  a  significant  impact  on  the  burning  rate. 


6.1  A  Fractal  Hypothesis 

In  our  earlier  discussion  of  the  planar  model  we  identified  the  well-known  results  that  the  fast- 
chemistry  limit  gives  rise  to  a  d2  —  t  law,  whereas  the  slow-chemistry  limit  gives  rise  to  a  d  —  t  law. 
Experiments  on  sub-micron  burning,  however,  give  laws  that  differ  subtantially  from  these.  Thus 
Parr  et  al,  [17]  report  a  d03  —  t  law  for  combustion  in  steam;  and  Allen  et  al,  [18]  report  a  d035  —  t 
law  for  air,  a  d0,24  —  t  law  for  a  CO2/H2O/N2  mixture. 

Parr  et  al,  [17],  discuss  agglomeration  problems  for  their  burner  experiments,  and  we  understand 
from  anonymous  sources  that  some  believe  that  agglomeraton  could  be  responsible  for  these  small-n 
laws.  But  one  would  think  that  the  agglomerations  would  need  to  be  fractal  in  nature  to  generate 


14 


Figure  3:  cj>  vs  a/b 


Hoop  Strains 


Figure  4:  Strains  vs.  a/b 


16 


Hoop  Stresses 


Figure  5:  Hoop  stresses  vs.  a/b 
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such  an  effect.  Allen  et  al,  [18],  examine  drop  burning  behind  a  reflected  shock  wave,  and  Allen 
(private  communication)  concludes,  using  conservative  collision  estimates,  that  agglomeration  is 
negligible  in  that  context.  And  so  here  we  propose  a  different  mechanism. 

It  is  reasonable  to  suppose  that  the  large  hoop  stresses  generated  as  the  aluminum/ alumina 
front  moves  inwards  will  lead  to  extensive  cracking.  It  is  well  known  that  large  cracks  in  metals 
can  be  fractal  in  nature,  and  one  might  speculate  that  this  is  also  true  for  the  cracks  in  the  small 
particles  of  interest  to  us.  Should  that  be  the  case,  the  front  will  be  fractal  in  nature. 

Suppose  we  are  in  the  slow  chemistry  regime  and  the  radius  of  the  aluminum  core  is  s  which 
changes  with  time  according  to 


ft=-K({Z},D,Trf)  (51) 

where  {H}  is  the  complete  set  of  parameters  that  characterize  the  ambient  conditions  and  the 
properties  of  the  aluminum  and  the  alumina,  and  Trf  is  the  temperature  at  the  reaction  front.  (We 
used  a  special  case  for  K  in  section  3  in  which  the  ambient  state  controls  c0,  which  in  turn  controls 
crf ).  We  assume  that,  overall,  the  dependence  of  K  on  s  (via  Trf  for  example)  is  small. 

Equation  (51)  can  be  written  as 


d_ 

dt 


=  -4tts2K 


(52) 


which  expresses  the  obvious  fact  that  the  rate  of  change  of  the  aluminum  volume  is  proportional 
to  the  aluminum  surface  area.  This  is  valid  for  a  sphere  with  a  smooth  surface  area.  But  if  the 
alumina  is  cracked,  aluminum  will  be  forced  into  the  cracks  by  pressure  gradients  and,  perhaps, 
surface  tension  effects  and,  since  the  cracks  would  be  thin,  this  would  substantially  increase  the 
aluminum/ alumina  contact  area.  If  s  is  now  redefined  to  be  the  mean  volume  radius 


where  V  is  the  total  aluminum  volume,  equation  (45)  can  be  rewritten  as 


d_ 

dt 


-4tt  s2+uK 


(53) 


for  some  parameter  v  which  accounts  for  the  enhanced  area.  This  choice  is  motivated  by  the  concept 
of  a  fractal  area  (in  3D  space)  of  dimension  2  +  v,  0  <  v  <  1  in  which  as  v  goes  to  zero  the  area 
dimension  approaches  that  of  a  smooth  surface,  and  as  v  goes  to  1  the  area  is  so  rough,  on  a  wide 
range  of  scales,  that  its  dimension  approaches  that  of  a  volume. 

Equation  (53)  implies  a  d1_l/  —  t  law  so  that,  for  example,  v  =  0.65  leads  to  a  d0'35  —  t  law. 


7  Conclusions 

In  this  paper  we  have  examined  a  simple  one-dimensional  model  problem  of  aluminum  combustion 
within  the  framework  of  the  shrinking  core  model;  and  the  spherical  problem.  The  conclusions 
of  the  first  problem  are  summarized  in  Section  5,  and  we  note  here  that  the  one-dimensional 
framework  permits  a  full  discussion  of  the  role  of  v,  the  speed  generated  within  the  oxide  by  the 
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O  and  A1  reaction,  together  with  unsteady  effects.  Such  matters  can  not  be  considered  for  the 
spherical  problem,  because  v  can  not  be  discussed  without  accounting  for  the  large  strain  fields 
that  it  generates.  However,  using  the  linear  equations  of  elasticity  we  can  estimate  these  strains  if 
we  assign  a  speed  to  the  reacting  front.  We  conclude  that  internal  strains  generated  thereby  are 
sufficient  to  crack  the  oxide  interior.  If  we  then  assume  that  the  crack  surfaces  are  fractal  in  nature, 
we  are  led  to  burn  laws  dn  —  t  for  which  n  can  be  significantly  smaller  than  1. 
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